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ABSTRACT 

We study the effect of spiral arm passages on the evolution of star clusters on planar 
and circular orbits around the centres of galaxies. Individual passages with different 
\Jq ■ relative velocity (Thrift ) and arm width are studied using iV-body simulations. When 

| the ratio of the time it takes the cluster to cross the density wave to the crossing 

time of stars in the cluster is much smaller than one, the energy gain of stars can be 
predicted accurately in the impulsive approximation. When this ratio is much larger 
than one, the cluster is heated adiabatically and the net effect of heating is largely 
C • ■ damped. For a given duration of the perturbation, this ratio is smaller for stars in the 

outer parts of the cluster compared to stars in the inner part. The cluster energy gain 
due to perturbations of various duration as obtained from our ./V-body simulations 
O is in good agreement with theoretical predictions taking into account the effect of 

adiabatic damping. Perturbations by the broad stellar component of the spiral arms 
on a cluster are in the adiabatic regime and, therefore, hardly contribute to the energy 
gain and mass loss of the cluster. We consider the effect of crossings through the high 
^ | density shocked gas in the spiral arms, which result in a more impulsive compression 

of the cluster. The energy injected during each spiral arm passage can exceed the 
total binding energy of a cluster, but, since most of the energy goes in high velocity 
k>( \ escapers from the cluster halo, only relatively little mass is lost. A single perturbation 

■ that injects the same amount of energy as the initial (negative) cluster energy causes at 

most ~40% of the stars to escape. We find that a perturbation that delivers ~10 times 
the initial cluster energy is needed to completely unbind a cluster in a single passage. 
The net effect of spiral arm perturbations on the evolution of low mass (< 100 M Q ) star 
clusters is more profound than on more massive clusters (> 1000 M Q ). This is due to 
the fact that the crossing time of stars in the latter is shorter, causing a larger fraction 
of stars to be in the adiabatic regime. The time scale of disruption by subsequent 
spiral arm perturbations depends strongly on position with respect to the radius of 
corotation (-Rcr), where Vdrift = 0. The time between successive encounters scales 
with Vd r ift as idrift <x Vdrift -1 and the energy gain per passage scales as AE oc V^ rift . 
Exactly at i?cR passages do not occur, so the time scale of disruption is infinite. The 
time scale of disruption is shortest at ~0.8-0.9 Rcr, since there Vdrift is low. This 
location can be applicable to the solar neighbourhood. In addition, the four-armed 
spiral pattern of the Milky Way makes spiral arms contribute more to the disruption 
of clusters than in a similar but two-armed galaxy. Still, the disruption time due to 
spiral arm perturbations there is about an order of magnitude higher than what is 
observed for the solar neighbourhood, making spiral arm perturbations a moderate 
contributor to the dissolution of open clusters. 

Key words: methods: iV-body simulations - galaxies: star clusters, spiral - Galaxy: 
open clusters and associations: general 
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1 INTRODUCTION 

The distinct difference between open and giobular ciusters 
has vanished since the discovery of yo ung massive ctusters in 
merging and inte racting gaiaxies (e.g. iHoltzman et at.ll 19921 : 
IWhitmore et al.l 1 1999ft . There is, however, stiil an evident 
difference in evolution. Star clusters formed in discs, which 
in our Galaxy are referred to as open clusters, experience ex- 
ternal perturbations by giant molecular clouds (GMCs) and 
by spiral arms and other disc density perturbations. These 
are not present in the halo of a galaxy, where most of the 
globular clusters reside. To understand how cluster popula- 
tions, such as the open clusters in the solar neighbourhood 
(|Kharchenko et al.| [2005) and such as the o nes found in spi- 
ral galaxies like M 51~ i|Bastian et al1l2005h and NGC 6946 
l|Larsen et al.ll200lh . evolve, it is important to understand 
the effect of these external perturbations. 

The number of Galactic open clusters as a function of 
age s h ows a lack of old open clusters, first pointed out by 
lOortJ (|1958T ). This lack can partially be explained by the 
rapid fading of clusters with age due to stellar evolution, 
which makes it harder to observe them at older ages. Still, 
fading can not explain the difference between the observed 
and the expected number of old (> 1 Gyr) open clusters, im- 
plying that a significant fraction must have been destroyed 
l|Wielenlll97ll ; lLamers et alj|2005h . In addition to two-body 
relaxation, clusters in tidal fields dissolve by the combined 
effect of: A.) Mass loss due to stellar evolution, reducing 
the mass and hence the binding energy of the cluster; B.) 
the tidal field of the host galaxy, imposing a tidal boundary, 
increasing the escape rate of stars; C.) bulge/disc shocks, 
pushing stars over the tidal boundary and D.) additional 
perturbations induced by irregularities in the galaxy, such 
as GMCs and spiral arms. The importance of the first three 
effects has been studied in detail by many people (for ex- 



ample IChernoff fc Weinberj|l990l:lGnedin fc Ostriker|[l997 



Gnedin, Lcc & Ostrikcr 1999; Takahas hi fc Portegies Zwartl 
200d : iBaumgardt fc Makind 120031 '). using different tech- 



niques (e.g. Fokker-Planck calculations or TV-body simula- 
tions). 

These studies were mainly aimed at understand- 
ing the evolution of globular clusters residing in the 
Galactic halo. The observed short disruption time (~ 

100 Myr) of disc cluste rs in the grand-design spiral 

galaxy M51 l|Gieles et all 12005ft is an order of magnitude 
lower than expected from the tidal fi eld of the galaxy 
l|Lamers. Gieles fc Portegies Zwartl I2005T ). Even for clusters 
in the solar neighbourhood, TV-body simulations predict dis- 
ruption times five times longer than observed. These TV- 
body simulations, however, ignore the presence of s piral 
arms (for example [T akahash ifc Portegies Zwartl |2000| and 



IBaumgardt fc Makind l2003h . In this paper, we study the 
contribution of spiral density waves to the cluster dissolu- 
tion time. 

Spiral arms are believed to be density waves rotating 
around the galaxy centre with an angular pattern speed (fi p ) 
whi ch is independent of the distance R to the galactic centre 
fsee lAthanassoulall984l for a review). The radius at which Q p 
is equal to the angular velocity in the disc (f2( r ) = Vdisc/R, 
where Vdisc the circular velocity in the disc) is called the 
corotation radius (Ren)- As gas travels through such a den- 
sity wave, it gets shocked and compressed to five to ten times 



higher densities !|Robertslll9"69T) . This is observable as nar- 
row sharp dust lanes in optical images of (nearly) face-on 
spiral galaxies which are far from edge-on. 

A cluster on a circular orbit around the centre of the 
galaxy will, in the inertial frame, have a velocity equal to 
the circular velocity Vdisc- In the reference frame of the spiral 
arm, its velocity (hereafter drift velocity, Vdrift) is equal to 



Vdrift = Vdisc — fipR 

= Viisc [1 - R/Rcn] ■ 



(1) 



The absolute value of Vdrift decreases going from the 
galaxy centre to -Rcr, where it is zero, and increases again 
beyond that radius. The time it takes for a cluster to travel 
from one arm to the next is defined as the drift time: 
tdrift(R) = 2nR/(m\VdTi!t(R)\), where m is the number of 
spiral arms in the galaxy. Thus, near to Rcr, a cluster ex- 
periences few passages but of long duration, while at con- 
siderably smaller, or considerably larger radii it undergoes 
many short-lasting spiral arm passages with high velocity. 
As it moves from the low density inter-arm region to the 
high density in the arms, the cluster gets compressed due 
to tidal forces. These accelerate the cluster stars, some of 
which may reach the escape velocity and become unbound. 

This is comparable to what happens to a 
globular cluster that crosses the Galactic disc 
l|Ostriker. Spitzer fc Chevalieij 1 1972ft . There are, how- 
ever, important qualitative differences between spiral arm 
perturbations and disc shocks. In particular, globular 
clusters cross the disc with a velocity almost independent of 
R. On the contrary, the velocity of the cluster with respect 
to the spiral depends strongly on the distance to the centre. 
Similarly, the scale-height of a dis c is independent of R, 
provid ed the gal axy is not barred dvan der Kruit fc Searld 
Il981al lbl ll982al lbT: Ide Griis fc van der KruitJ 1 1996ft . For 
spiral arms, on the contrary, it is the angular width of 
the arm that is roughly independent of radius, so that the 
linear width of a spiral arm is proportional to the radius. 

The passage duration determines the nature of the per- 
turbation of the arm on the cluster. If the time to cross the 
arm (T cr ) is much longer than the crossing time of stars in 
the cluster (t CI ), the stars can adiabatically adjust to the 
density increase. When T CT is much shorter than t cr , the 
stars get an impulsive velocity increase due to tidal forces. 
Stars in the cluster core have a short t cr and, therefore, get 
adiabatically heated. The boundary between adiabatic and 
impulsive heating in the cluster is determined by the ratio 
tcr /Tcr- If this ratio is low, a higher fraction of stars is heated 
impulsively. 

In this paper we study the effect of spiral arm passages 
on clusters and investigate the importance of relative veloc- 
ity with respect to the arm. This paper is organised as fol- 
lows: In Section[2]the physical parameters of spiral arms are 
derived from observations. Simple analytical estimates of the 
effects of compressive perturbations are given in Section [3] 
They are confronted with results of TV-body simulations in 
Section[4] In Section[5]the results of the simulations are used 
to derive cluster mass loss and disruption times due to the 
spiral arm perturbations. The conclusions are presented in 
Section [6] 
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2 SPIRAL ARM PARAMETERS 
2.1 Stellar arms 

Disc clusters are found to have short disruption times both 
in our Galaxy and in M51. Therefore, we will hereafter 
consider two fiducial spiral galaxies, one Milky- Way-like - 
which we will call for brevity Milky Way, or MW - and 
the other M51-like - which we will call grand design spi- 
ral. We will adopt for both a constant rotational velocity of 
Vdisc = 220 kms" 1 , that is, we limit ourselves to the regions 
where their respective rotation curves are flat. 

Spiral patterns are presen t in both youn g and old pop- 
ulations, as alre ady noted by IZwickvl l|l955ft and later by 
ISchweizer] l|l976T ). The old Population II stars in the disc 
show a broad spiral pattern. In contrary, the young Popu- 
lation I stars, which form preferentially in the spiral arms, 
show more irregular spiral patterns. Estimates of the arm- 
interarm density contrast (r) of the stellar p o pulati on in 
th e disc of M51 were made by iRix fc Riekj (J1993) and 
by iGonzalez fc Graham! (|l996l ). using the A"-band light as 
tracer of stellar mass. Their results are in good agreement 
and show that the maximum value for P is around 2 — 3. 
Similar values have been found for other grand d esign spi- 
rals (e.g. NGC 4254, see lKranz. Slvz fc RixllioOll ). We ex- 
pect our own Galaxy to h ave smaller values of P. Indeed, 
iDrimmel fc Spergell (|200lh find a value of F ~ 1.3 from 
model comparisons to K-band measurements. 

M51 has a clear two-armed (m = 2) grand-design spi- 
ral, presumably due to the interaction with its compan- 
ion NGC 5195. The stru cture of our G alaxy is more com- 
plex. It has a bar (see iDehnenl |2002| for a review and 
references therein) and a spiral structure which could ro- 
tate with a p attern speed differe nt from that of the bar 
dTagger et all 1 19871; ISvgnet et al.1 1 19881 ). iDrimmell (j2000l ) 
and lDrimmel fc Spergell (|200lft modelled the COBE-DIRBE 
data and propose that the old stellar population has a two 
armed spiral, while the gas and youn g and stars form a four- 
arme d spiral (for the latter see also iGeorgelin fc Georgelinl 
1 19761 ). Studies in which all components (bar and spirals) 
have the same pattern speed fin d that the sola r neighbour - 
hood is well beyond corotation l|Dehnenll2OO0l ; lFuxll200lft . 
but models allowing for a different pattern speed for the 
bar and spiral components gi ve a better view of th e com- 
plex structure in our Galaxy. iBissantz et al.1 (|2003l ) calcu- 
late the gas response to a composite bar plus spiral model 
and find best agreement for the observations with a pattern 
speed of 60 kms" 1 kpc" 1 for the bar and 20 kms" 1 kpc" 1 
for th e spiral. Similar values are found by iMartos et al.l 
|2004l ). This places the solar neighbourhood well beyond the 
bar corotation, but just within the spiral corotation. With 
an adopted solar radius of Ro = 8.5 kpc and the adopted 
value for Vdisc = 220 kms" 1 the solar neighbourhood is at 
R ~ 0.8i?cR- Compared to the Milky Way, M51 has been 
modelled very little and no clea r-cut value has been given 
for its patte rn speed. Recently, Zimmer. Rand fc McGrawl 
(|2004 ) and lEgusa. Sofue fc Nakanishil (|2004l ) determined 
the pattern speed of four grand design two-armed spirals 
and they find values, between 30 and 40 kms -1 kpc -1 . 
We adopt Rck~ 6 kpc, resulting in a pattern speed of 
~ 37 kms" 1 kpc" 1 . 

To further quantify the effect of a spiral-arm passages 
on the evolution of star clusters, we also need to consider 



the width of the arms. To enable a quantitative discussion 
about spiral-arm width, we define A <?5>i/2 as the arm an- 
nular full width at half maximum. The typical values of 
A0!/ 2 measured from optical images of six spiral galax- 
ies by ISchweizel (Il976l). is almost in dependent of R and 
A cj> 1/2 ~ 20°. ISeigar fc James! (| 1998ft studied K-band im- 
ages of a sample of 45 face-on spiral galaxies and found for 
two-armed spiral galaxies typical values of A4>i/2 between 
20° and 40°. The se values are slightly higher than found by 
ISchweizerl l| 1976ft . probably because the A"-band is a bet- 
ter tracer of the old stellar populations, for which 
is larger. For A(/>i/2 = 30°, the value of T cr depends on 
Vdrift and R as T cr ~ 0.5-R/j Vdrift | - For the solar neighbour- 
hood T cr — 80 Myr. For our M51-like galaxy, the values for 
T cr at [0.3/0.5] x R C n are [14, 68] Myr. The values for T cr 
are all higher than the typical t CI of clusters, that is, a few 
Myrs. Therefore, the stars in the cluster respond adiabati- 
cally when the cluster crosses the arm and this crossing does 
not have the short lasting compressive effect as a galactic 
disc has on globular clusters. 

2.2 Gaseous arms 

The gaseous spiral pattern is different from the stellar one. 
Gas gets shocked in the potential well of the stellar arms, 
whe re it gets com pressed to up to ten times higher densi- 
ties l|Robertslll969ft . These high gas densities manifest them- 
selves in optical observations of spiral galaxies as thin dust 
lanes on the leading side of the arm within Ren- After 
the shock, the density decreases exponential-like to a value 
slightly lower than the mean gas density, as is illustrated 
in Fig. [T] Such strongly peaked high gas densities act as a 
compressive perturbation on the cluster (Section which 
can be compared to disc shocks of globular clusters. The disc 
crossing is referred to as a shock in that context because of 
the short duration of the perturbation. That is, the time to 
cross the disc is much shorter than t CT . Since we here also 
discuss shocks by gas in the arm we avoid confusions and 
refer to perturbations instead of shocks. In addition, in Sec- 
tion |3] we will show that, due to the variation of T cr with R, 
the perturbations are of highly varying duration. 

Besides the higher densities, also the typical values of 
A (/>i/2 of the ga seous spiral arm are much smaller than the 
one of the stars. iRobertsI (|l969l ) argues that the width can 
be calculated by the time scale of formation and evolution of 
massive stars, that is, a few tens of Myr, since these stars are 
visible in the optical, so they must have traveled out of the 
high density dusty environment where they for med. Studies 
of the mole cular spiral arms in M31 and M51 bv lNakai et al.l 
1 19941 ) and lNieten et al.1 l|2006l ). give values of the full- width 
at half maximum (W1/2) of 500 pc and 1000 pc, respectively. 
Since these values are close to the resolution of the beam, the 
sharp peak in density (Fig. [l]) is presumably not resolved. 
Hydrodynamical simulations of barred galaxie s revealed nar- 
rower density peaks of the gas (see Fig. 11 of I Athanassoulal 
1992), very close to what was predicted bv lRobertsI (|l969D 
and what is shown in Fig. [T] We also measured the width of 
the dust lane on an optical picture of M5 1 and the w i dth o f 
the HI on the 8 arcsec resolution map of lRots et ail (|l990l ) 
and find that they are both compatible with a width of 250 
pc. We, therefore, adopt Wi/ 2 = 250 pc for the gaseous arms. 

Measurements of the neutral hydrogen content in a sam- 
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Table 1. Adopted parameters for spiral arms in our fiducial grand 
design two-armed spiral galaxy and in the Milky Way 





Grand design 


Milky Way 


rn 


2 


4 


Vdisc (kms -1 ) 


220 


220 


Op (km s — 1 kpc — 1 ) 


36.7 


24.5 


Rcb. (kpc) 


6 


9 


A </>i/2 (stars) 


30° 


15° 


W"i/ 2 (gas) (pc) 


250 


250 


Tstars 


2-3 


< 2 


Tgas 


10 


10 


(Phi) (M q pc" 3 ) 


0.033 


0.033 



pie o f spiral galaxies in the Virgo cluster l)Cavatte et al.l 
1 19941 ) show that the surface density of HI gas (Shi) within 
R ~ 10 kpc in gas rich spirals is nearly constant with R 
and of the order of Shi — 10 Mqpc~ . A recent study of the 
Milky Way H cvinc, Blitz fc He ilcs 2006) gives similar v alues 
for Shi in the solar neighbourhood. iLevine et al.l (|2006h also 
find that T is almost independent of R. Their measurements 
of the arm-interarm surface density contrast, combined with 
the measured vertical compression factor in the spiral arms 
are consistent with V = 10 for the gas. 

From Shi and the vertical scale height of gas {H z ) we 
can derive the mean midplane density of HI ({pm)), since 
(pm) = Shi/2_ Hz, assuming an exp onential vertical den- 
sity distribution. Ivan der KruitJ |l98Sl) finds a vertical scale 
height of H z ~ 150 pc, resulting in (pm) — 0.033 M Q pc -3 . 
The peak density (po) is then 

po-r<pm), (2) 

where we assume F is the ratio of po to (pm). 

All parameters of spiral arms, used in our models, are 
summarised in Table [T] 

2.3 Potential-density relations for the arms 

The density of gas along the orbit of the cluster has a sharp 
rise at the point where the shock occurs and an exponential- 
like decrease after the shock (see our Fig. \T\ and Fig. 5 of 
I Roberts! 19691) . In reality, the density increase is bit smoother 
than shown in Fig. [T] It is hard, if not impossible, to derive 
a suitable potential that describes a realistic density func- 
tion. Therefore, we assume for simplicity a one-dimensional 
density perturbation of Gaussian form 

p(Y) = po exp(-Y 2 /# 2 ), (3) 

where H is the scale length and Y is the azimuthal distance. 
In Section [3] we show that for an impulsive perturbation the 
compressive tidal forces on the cluster scale with the integral 
of p(Y) with respect to Y, where Y is the direction of motion 
of the cluster. A Gaussian function with a integral equal to 
that of an exponential function has the same compressive ef- 
fect on the cluster. We emphasise that this is only true when 
for both perturbations hold that they are impulsive. In prac- 
tise, this implies that the majority of the surface should be 
contained within the same Y extent for both functions. That 
is, a very broad gaussian with the same surface as a very 




-2 2 4 6 

Y/W,/2 

Figure 1. Schematic illustration of the gas density along az- 
imuthal trajectories through a spiral arm, where the gas travels 
from left to right in this figur e . The grey shaded area is the gas 
density predicted by iRobertj jl969h . When the scale length of 
a Gaussian function (H) is times the scale length of the 

exponential function (Wi/2)^ t ne areas under the two curves are 
the same. 



peaked gaussian will have a different effect on the cluster, be- 
cause the broad gaussian will heat the cluster adiabatically, 
while the peaked gaussian will heat the cluster impulsively. 
For the exponential function p(Y) = po exp(— Y/W1/2), for 
Y > 0, and the Gaussian function of equation (O, the sur- 
faces are equal when H = W1/2/ 'y/ft- With Table Q] we thus 
find that H ~ 150 pc for the gaseous arms. 

By using Poisson's law and equation <(3j we can derive 
the related potential along the trajectory of the cluster 

4>(Y) = 2-nGp Q H [v^rYerf (Y/H) + H exp(-F 2 /^ 2 )] > (4) 

where G is the gravitational constant. The acceleration 
(g(Y)) that is felt by the cluster due to spiral arm passage 
is equal to — dcj>/dY and with equation Q it follows that 

g(Y) = -2n 3/2 Gp Hetf(Y/H). (5) 

The Gaussian density form (equation |3} implies a constant 
acceleration far from the centre of the density wave. This is 
not physical for spiral arms, since the attractive force due 
to the arm should decrease with distance from it and should 
reach zero in between two arms. Therefore, we only consider 
the effect of the density wave in the vicinity of its centre 
(-2H < Y < +2H), that is, where g(Y) varies and hence 
tidal forces are at work causing a compression of the cluster 
(Section 0). 

The values of A $1/2 and W1/2 are always measured in 
the tangential direction, that is, along the orbit of the clus- 
ter. The width perpendicular to the spiral arm is therefore 
smaller and depends on the pitch angle. Given the uncer- 
tainties in the arm width measurements and since the typi- 
cal radius of star clusters (~ 3 — 4 pc) is much smaller than 
the W1/2 of the spiral arm, it does not matter that clusters 
cross the arm with a certain angle. In order to avoid having 
the pitch angle as an extra parameter, we study perpendicu- 
lar passages through a density wave. The effect of the pitch 
angle is taken into account in the choice of the azimuthal 
scale length H. 
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3 SIMPLE ANALYTICAL ESTIMATES FOR 
ONE-DIMENSIONAL TIDAL 
PERTURBATIONS 

3.1 The impulsive approximation 

The energy gain of a globular cluster cro ssing the Galactic 
disc was derived bv lOstriker et afl (ll972T ) using the impulse 
approximation. They assumed that the stars do not move 
during the passage of the density wave, that is, that t CI is 
much longer than T cr . 

Assume a one-dimensional acceleration along the Y-axis 
(g(Y)), a cluster with its centre at Y c and an individual star 
at a distance y = Y — Y c from the cluster centre. The tidal 
acceleration of a star due to the density wave is then 

dv v 



: g(Y) - g(Y c ) ~ y 



(6) 



d<? y dg 

At ' »^'-»dy = vdt' 

where we have substituted dY = Vdt, with V the relative 
velocity between the cluster and the density wave and where 
g is expanded around the cluster centre. This is sufficiently 
accurate as long as the cluster is much smaller than the 
width of the density wave. Note that the tidal acceleration 
scales with the density variation (p(Y)) since g — — d(f>/dY 
and p(Y) oc d 2 4>/dY 2 through Poisson's law. Therefore, the 
tidal forces scale with — p(Y) (equation [6]) and are always 
directed inwards for a density increase. Com bined with the 
impulsive assumption, lOstriker et al.l (|l972h introduced the 
name compressive gravitational shock. As mentioned in Sec- 
tion [2l the term shock was introduced due to the short du- 
ration of the perturbation. This does not necessarily have to 
be true in our case, since for cluster close to Rcb. T ct can be 
long. 

Integrating equation ([6| then yields an expression for 
the velocity increment of a star of the form 



y\9(Yc) \_y ni 



where Y m is the point where g reaches its maximum g m . If g 
is an odd function and does not change much for Y < — Y m 
and for Y > +Y m , the total energy gain per unit cluster 
mass after a tidal perturbation is 



AE iu 



2»„ 



3 V 2 



1„2 



(8) 



and r 2 is the mean 



where we have substituted y 2 
square position of stars in the cluster. From equation ([5]) it 
follows that gr m is 



g m = 2n 3/2 Gp H 



(9) 



As mentioned in Section 12.31 the constant acceleration at 
large distances from the Gaussian density is not physical 
in the case of spiral arms. However, we are interested in 
the density change and the related change in g from — <? m to 
+gm close to the centre of the density wave is approximately 
correct. 



3.2 Validity of the impulsive approximation 

3.2.1 Constant velocity assumption 



One of the assumptions made bv lOstriker et ail (| 19721) and 
later in more thorough studies ( see for example IWeinberd 
Il994d ; iMurali fc Weinberg! Il997i ; iGnedin fc Ostrikerl 1 19971 




-0.5 0.0 
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Figure 2. The variation of 1/V with g(Y) (equation |6j for a 
constant velocity equal to Vo and V ma x (dotted and dashed, re- 
spectively) and for the true velocity V(g) (full line), derived from 
the <f>(Y) (equation [1J. 



1999) on the tidal perturbation due to the Galactic disc 
on globular clusters is that the velocity, V(Y), remains con- 
stant during the crossing. This is probably not such a bad 
assumption, since globular clusters cross the disc with high 
initial relative velocity (~ 200 kms -1 ). For spiral arm cross- 
ings, however, the relative velocity is between zero at Rcr 
and almost il^Usc close to the galaxy centre and at large 
distances from it. For low | Vdrift | there is large relative in- 
crease of velocity, making the assumption of a constant ve- 
locity invalid (equation |S} and|S}. The solution to this can 
be obtained by expressing V as a function of g and inte- 
grating equation ©. This, however, is even for this simple 
functional form quite hard. 

Alternatively, we can estimate the change in V(Y) from 
4>{Y) (equation Q. The variation of V(Y) depends on <f>(Y) 



V(Y) = V / Y O 2 + 2[^ O -0(Y)], 



(10) 



wherex cj>o and Vo are a reference potential and velocity, 
respectively. The variation of g(Y) follows from the spatial 
derivative of (f>(Y). We need to integrate 1/V with respect 
to g (equation [6]), to get the total Av y of an individual star. 
In Fig.[2]we show the variation of 1/V with g for a constant 
V (dotted) line and one that considers an acceleration due 
to the potential (full line). The shaded surface is the result 
of the integration of 1/V with respect to g. This surface is 
almost equal to ~ 2g IU /V lnax . 

Using equation Q, we can find an expression for A</> 
for Y » H of the form 



A<t> = 4>{Y ) ~ 4>(0) ~ 2ttG Po H [V^Yq - H] 



(11) 



where Yq is the starting position of the cluster, where it has 
velocity Vo with respect to the density wave. Equation (|I1[) . 
combined with equation QD) , can be used to find the value 
of Vmax. This suggests that it is better to use Vmax instead of 
Vo in equation (JSj to calculate the energy gain. In Section |4] 
we validate this with TV-body simulations. 
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3.2.2 The effect of adiabatic damping 

The value of t CI depends on the distance of the stars to 
the cluster centre (r). In the core of the cluster t C r is much 
shorter than close to the tidal radius (rt). For stars with 
short period orbits, the effect of the shock is largely damped 
due to adiabatic invariances. Therefore, there is almost no 
net energy gain in the core. On the other hand the stars in 
the outer region will largely be heated impulsively. To define 
the transition between the impulsive a nd adiabatic regions, 
we can define an adiabatic parameter (|Spitzerlll987h 

(12) 

where u) is the angular velocity of stars inside the clusters. 
This is defined as u)(r) = <x(r)/r, with a(r) the velocity 
dispersion of stars at position r. (Not e that the orig inal def- 
inition of the adiabatic parameter bv lSpitzerlll987l contains 
an additional factor 2. We choose to define x as the ratio of 
the time it takes the cluster to cross a distance H to the time 
it takes a sta r to cross a distance r in side the cluster, e.g. 
T a /tct, as did lGnedin fc Ostrikerl l|l999l U When x « 1 the 
term shock is justified, while for x >> 1 the perturbation is 
larg ely adia b atic. 

ISpitzerl (|l95St ) gives an estimate of the conservation of 
adiabatic invariants in the harmonic potential approxima- 
tion. This assumes all stars are initially at the same distance 
r from the cluster centre and have the same oscillation fre- 
quency co. The energy gain for each star can then be written 
as a function of the result of the impulsive approximation 
fSection l3.ll) and an adiabatic cor rection factor A(x) (equa- 
tion 5 of Kundic & Ostrikedll995h 



AE = AE imp A(x) 



(13) 



where A_E"i mp is given by equation (JS| . The correctio n factor 
A(x) in the harmonic approximation is (for example ISpitzerl 
1 1981 Eqs. 5-28) 

2\ 



A s (x) = exp(— 2x 



(14) 



We refer to A s (x) as the Spitzer correction. (The original 
definition of A s (x) by Spitzer has a factor 0.5 instead of 2, 
bec ause of the differen t definition of x in equation (|12p .) 

IWeinberd l|l994aH blloh showed that this exponential de- 
crease of the energy with x underestimates the heating ef- 
fect. This is due to the fact that the basic assumption of 
the harmonic potential approximation is not valid when the 
system has more than one degree of freedom, so that small 
perturbations can still grow. When the cluster is represented 
as a multidimensional system of nonlinear oscillators, some 
perturbation frequencies become commensurable with the 
oscillator frequencies of stars. This results in a correction 
factor A(x) that is not exponentially small for large x but, 
instead, has a pow er-law dependence on x. The simplest 
form, as shown bv lGnedin fc Ostrikerl l|l99l ). can be writ- 
ten as 



A w (x) = (l + x 2 ) 



-3/2 



(15) 



We refer to A w (x) as the Weinberg correction. 

Besides this shift in energy (AE), there is also a 
quadratic term that affects the dispersion of the energy spec- 
trum of the clust er. The effect of shock-ind uced relaxation 
was mentioned bv lSpitzer fc Chevalier! (Il973f) and later stud- 
ied in more detail by Kundic fc Ostriker (|l995h . It is beyond 



the scope of this study to include all this detailed physics 
and we refer the reader to the aforementioned papers for 
details. 

In Section 0] we will compare results from iV-body sim- 
ulations to both adiabatic correction factors (Eqs. [14] and 
[15]). 



4 A-BODY SIMULATIONS 

In this Section, we confront the simple analytical estimates 
given in Section [3] with a series of iV-body simulations of 
clusters that cross a density wave of Gaussian form (equa- 
tion [3]). The potential and acceleration are derived from the 
one-dimensional Gaussian density profile, as described in 
Section 12.31 



4.1 Set-up of the simulations 

4-1.1 Description of the code 

The A-body calculations were carried out by the kira in- 
tegrat or, which is part of the Starlab software environ - 
ment l|McMillan fc HutJll99d : iPortegies Zwart et al.ll200lh . 
Kira uses a fourth-order Hermite scheme and includes spe- 
cial treatments of close two-body and multiple encounters 
of arb itrary complexity. T he special purpose GRAPE-6 sys- 
tems ( Making et ail 120031 ) of the Observatoire de Marseille 
and of the University of Amsterdam were used to accelerate 
the calculation of gravitational forces between stars. 



4-1.2 Units and scaling 

The cluster energy per unit mass (-Eo) is defined as Eo = 
7/GM c /(2rh), with r\ ~ 0.4, M c is the mass of the cluster 
and rh is its half-mass radius, depending on the cluster 
model. All clusters are scaled to A-bo dy units, such that 
G = M c = 1 and E = -1/4, following iHeggie fc Mathieul 
( 1986). The virial radius (r v ) is the unit of length and follows 
from the scaling of the energy, since r v = GM C /(2\U\) = 1, 
where U is the potential energy per unit of cluster mass. 
We assume virial equilibrium at the start of the simulation, 
which implies U — 2Eq and, therefore, r v (t = 0) = 1. 



4-1.3 Cluster parameters 

For th e density distribution of the cluster we assume a lKind 
<l 19661) profile, which fits the radial luminosity profile of our 
Galactic clusters. Open clusters in our Galaxy have a lower 
concentration index than globular clusters. Defining the con- 
centration index as c = log(rt/r e ), where rt is the tidal ra- 
dius and r c is the core radius of the cluster, we find that typ- 
ical values f or open clusters are in the range 0.5 < c < 0.9 
(|King|ll962h or even lower l|Binnev fc Tremaind [l987h . A 
concentration index c = 0.7 corresponds to a dimension- 
less central potential depth of Wo = 3. The a verage concen - 
tration index of globular clusters is c ~ 1.5 (|Harris|[l996l '). 
corresponding to Wo = 7. Here, we adopt a dimensionless 
central potential depth of Wo = 5. This corresponds to a 
concentration index of c = log(r t /r c ) ~ 1.03. For a Wo ~ 5 
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Figure 3. Variation of fractional energy gain (AE/\Eq\), shown 
as a full line and of the kinetic and potential component shown 
as dotted and dashed lines, respectively. The cluster has N = 
65 536 particles and started in virial equilibrium. Parameters for 
the density wave are described in the text f Section 14. 1.41 . 

cluster it follows that rj — 0.41. In iV-body units, the cor- 
responding radii are r\ = 0.81 and r 2 — 1.36, respectively 
(Section ELS) . 

4-1-4 Code testing 

To test our code, we run a single perturbation with H = 
100 and Vo = 200 and po chosen such, that AE/\Eo\ — 
0.05 (Eqs. [8] and [9}. In this example, the cluster consists of 
N — 65 536 equal mass particles. It is initially positioned at 
Yo — —2H and its velocity is directed towards the density 
wave. When the cluster is at +2H, there are almost no tidal 
forces anymore, but there is still an acceleration towards 
the density wave (equation [SJ. To prevent a second crossing, 
we turn the external tidal field off. The cluster is evolved 
for an additional 10t C r after the perturbation. In Fig. [3] we 
show the variation of the energy gairQ (AE) as a full line 
and of the internal kinetic energy (AT) and potential energy 
(At/) as dotted and dashed lines, respectively. Due to the 
compressive nature of the tidal forces, all stars gain kinetic 
energy in the —y direction, resulting in an increasing T. 
Consequently, the stars move deeper in the potential of the 
cluster, resulting in a decreasing U. The cluster is in virial 
equilibrium before the perturbation, that is, the virial ratio 
is equal to one (Q — —2T/U — 1). The value of Q increases 
due to the perturbation. After the perturbation the cluster 
revirialises, due to relaxation, reducing Q initially to Q < 1, 
around t ~ 3t CI . After a few oscillations the virial ratio is 
close to one again, leaving a cluster with AT = — AE and 
AU = +2AE. 

4.2 Parameters of the runs 

We choose to use the same clusters as before, but now with 
iV=8k. The value of uj varies strongly within the cluster, 

1 When we discuss cluster energy, we refer to the sum of the 
internal potential and internal kinetic energy, that is, where the 
energy of the centre of mass motion and the contribution of the 
external tidal field have been subtracted. 




0.01 0.10 1.00 10.00 

r 



Figure 4. Frequency w (full line, scale on the left) and the mass 
contribution Airr 2 p(r) as a function of r (dashed line, scale on 
rigth). The value of equation H16I I is indicated with an arrow. 

from > 100 deep in the core down to ~ 0.01 at rt- We, 
therefore, calculate the mass weighted mean value of u), i.e. 
uj, by numerically solving 

47r r Tt 

Q = — J r 2 p(r)oj(r)dr. (16) 

For a King profile with Wo = 5, u> = 0.68. This is about 
a factor of two higher than l/t CI , with t CI = 2\/2 at r v (in 
iV-body units). The radius at which oj(r) — 0.68 is slightly 
inside the half-mass radius: r(uj = Q) = 0.72rh. Fig. f4] illus- 
trates the variation of uj(r) (dashed line) and on 4-7rr 2 p(r) 
(full line) for a Wo = 5 cluster. 

The value of po was chosen such that the predicted 
ABimp (equation[8| is always 1/8 of \Eo\, which is +1/32 (in 
iV-body units) (Section l4.1.2[) . This is possible since po is a 
free parameter in Eqs. [5] and [5] after Vo and H are defined 
by our grid. 

The cluster is positioned initially at a distance Vo = 
—2H with a velocity Vo in the Y-direction. 

4.3 Testing the validity of the impulsive 
approximat ion 

We ran two series of 12 simulations, both for the same se- 
quence of x values (equation [T^J , ranging from 0.003 to 15. 
This was achieved, in the first series by fixing Vo = 50 and 
varying H and in the second series by fixing H = 50 and 
varying Vo accordingly. Simple analytical estimates relying 
on the impulsive approximation (equation [Sj give that the 
energy gain AE scales as 1 /Vo 2 and is thus independent of 
x. If, however, we use in equation ((SJ 1/V^ ax instead, we 
find that AE decreases with increases x. 

In Fig. [5] we show AB/|i?o| as a function of x, defined 
in terms of Vo- The result for constant Vo and H are shown 
in the left and right panel, respectively. The square sym- 
bols are the results of the iV-body simulations. The dashed 
line shows the prediction of A£j m p (equation [SJ based on 
a constant velocity. Excellent agreement between the sim- 
ulations and the analytical predictions using the impulsive 
approximation is found for x < 0.03, for both series. With 
the full line we show the predicted energy gain when we 
use V m ax in equation ((8}, which we derived from Eqs. [10] 
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Figure 5. Fractional energy gain (AE/\Eq\) as a function of the adiabatic parameter x (equation I12t resulting from the iV-body 
simulations described in Section 14,21 (squares) , the impulsive approximation (equation [8} based on Vb (dashed line) and on V max (full 
line). In the left panel the value for Vb was fixed and H was varied to vary x. In the right panel it was done the other way around. The 
two series give the result. 



and 1111 Indeed the energy gain decreases for larger x, since 
Knax/Vb increases for perturbations with longer duration. 
This prediction of Ai?i mp (V max ) explains the results of the 
simulations very nicely for x < 0.5. For x > 0.5 the energy 
gain in the simulations is lower than predicted by the impul- 
sive approximations. This is caused by adiabatic invariances 
in the cluster core (Section 13,2.20 . The two sequence give 
nearly identical results, which suggests that it is the ratio 
H/Vo, or the duration of the perturbation that controls the 
energy gain together with po- In Section ^. 4l we compare the 
results of the simulations to the theoretical predictions that 
include adiabatic damping. 



4.4 Adiabatic invariances 

To quantify deviations from the impulsive approximation, 
that is, the effect of damping by adiabatic invariances, we 
show in Fig.[U]the resulting AE /\Eo\ of the simulations from 
Section |4]3j but now using Vrnax in the definition of x (equa- 
tion HU) and A_E imp (equation [HJ). The values of AE/\E \ 
from the simulations are multiplied by (V^nax/Vb) 2 . If all sim- 
ulations were in the impulse regime, the energy gain would 
be constant (AE/\Eo\ = 1/2), which is shown as a dashed 
line. (Note that the dashed line in Fig. |S] is the same as 
the full line in Fig. [5}. The theoretical predictions with adi- 
abatic correction are shown for the result of Spitzer (equa- 
tion !14p and Weinberg (equation ll5p as dotted and full lines, 
respectively. As was already clear from Fig. [5] the simula- 
tions can be well explained by the impulsive approximation 
for x < 0.5. The decrease of AE/\ E \ for x ^ 0.5 is well 
explained by the predictions made bv lWeinberel |l994al lbHcl). 
that is, equation (|15p (full line). 

From Table [1] we find that the values of x for a cluster 
with O = 1 Myr and stellar arms at R — [1,3, 5] kpc in a 
grand design spiral are ~ [3, 10, 50], respectively. Combined 
with the results of Figs.[5]and[6]this support our assumption 



that the stellar arms do not contribute to the disruption of 
clusters (Section [2}. 

For the gaseous arms close to Rcr the values of | Vdrift | 
are low, and x > 1 (Section 0). It is, therefore, important 
to include a correct description of adiabatic damping to un- 
derstand the effect of spiral arm perturbations on clusters. 
Including the simple correction as derived by Spitzer (equa- 
tion I14[) underestimates the energy gain of perturbations 
with low | Vdrift |- The impulsive approximation (equation [Sj 
overestimates the energy gain by orders of magnitude. 



5 CLUSTER MASS LOSS AND DISRUPTION 
BY SPIRAL ARM PERTURBATIONS 

5.1 Energy gain vs. mass loss 

The energy gained in a compressive perturbation is absorbed 
mainly by the stars in the outer parts, since for these stars 
r is large (Eqs. [6] [7]). This implies that the fractional en- 
ergy gain (A£V|_Bo|) is not necessarily equal to the frac- 
tional mass losfO {AM /Mo). In fact, AE/\E \ is higher than 
AM /Mo, since stars escape with velocities higher than the 
escape velocity. Since this energy is taken by escaping stars, 
it can not be used to unbind more stars. This w as also found 
for en counters between star clusters and GMCs (|Gieles et alj 
120061) . 

We simulate perturbations with A_E/|_Bo| between 10 -3 
and 100 for clusters with N = 2048. The values of H and Vb 
were fixed at 10 and 100, respectively. The value of po was 
varied to achieve the desired energy gain. The simulations 
are continued without a tidal field for five crossing times 

2 We define the mass loss, AM, to be positive. Bound is defined 
as having a velocity lower than the escape velocity, which is based 
on the cluster potential only, that is, corrected for the presence 
of the external potential. 
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after the perturbations. The final number of unbound stars 
is then compared to A_B/|_Bo| from the simulations. In Fig. [7] 
we show AM /Mo as a function of AE/\Eo\ following from 
the simulations (circles). The dotted line shows a one-to-one 
relation. For AE/\E \ < 1 the relation between AM /Mo 
and AE/\Eo\ is almost linear, but with AM/Mo a factor 
of 7.5 lower than AE/\E \. For AE/\E \ > 1 the relation 
flattens to AM /Mo = 1, that is, to a value corresponding to 
a completely unbound cluster. This relation can be described 
well by a simple function of the form 



AM/Mq = 1 - exp(-fAE/\E \), 



(17) 



with / ~ 1/7.5 (dashed line). The reason for the deviation of 
a one-to-one relation at AE/\Eq\ ~ 1 is that the assumption 
of the impulsive approximation is not valid anymore for such 
strong perturbations. The stars in the centre of the cluster 
are heated adiabatically (Fig. [4} and remain bound after the 
perturbation. A value of A_E/|_Eo| of almost ~10 is necessary 
to completely unbind the cluster with a single passage. 

Stars that do not get unbound can get in higher en- 
ergy orbits, where in the presence of a tidal field they can 
be beyond r t . In these simulations we ignore the Galactic 
tidal field, but we can follow the number of stars that is 
pushed over rt. For this, we adopt a physically more rele- 
vant definition of AM/Mo that considers all unbound stars 
as well as stars that are still bound to the cluster, but have 
positions larger than r t . The relation between this AM/Mo 
and A_E/|£"o| is shown as squares in Fig. [7] These values 
are higher, but still a factor of 2.5 lower than AE/\Eo\, 
and the results can be approximated by equation (|17[1 , with 
/ ~ 1/2.5 (full line). To translate the energy gains predicted 
from the analytical estimates to mass loss of the cluster we 
adopt equation (|17[) with / = 0.4. With this we estimate 
the mass loss and resulting disruption time of clusters due 
to spiral arm perturbations in Section f5. 21 



5.2 The cluster disruption time 

5.2.1 A toy model for spiral arms in the disc 

The strength of a spiral arm perturbation on a cluster de- 
pends on Vdrift, H and po. In Sections 12.21 and 12.31 and Ta- 
ble [1] we adopted a constant arm-interarm density contrast 
F = 10 and scale length of H — 150 pc for gaseous arms, 
leaving Vdrift as main R dependent variable. 

From equation ([2]) and (pm) from Section [2.21 we find 
po = 0.33 Mq pc -3 . The value of Vmax in the arm is found 
from Vdrift and by considering an acceleration by the arm 
from Yo = —2H to Y = (equation I10p . With these param- 
eters and Eqs. [13] and [15] we can now calculate the energy 
gain for clusters that cross spiral arms at different R. 

With the relation between AM/Mo and AE/\E \ from 
Section ["5.1l that is, equation (|17|l with the factor / = 0.4, we 
can estimate the mass loss of the cluster at different R. We 
assume here that AM/Mo is independent of N, which can 
be justified since the derivations of AE/\Eo\ (equations. [8] 
and I13[) are independent of N as well. They do, however, 
depend on the distribution of stars in the clusters, that is, on 
the cluster density profile or the concentration index. There 
might be a small N dependence in the boundary between 
the impulsive and the adiabatic regime, since clusters with 
low N are able to adjust faster to a changing potential. 
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Figure 7. Fractional mass loss defined as the fractional loss of 
bound stars (AM/Mq) vs. fractional energy gain (AE /\Eq\) from 
the simulations (circles). The dashed line is a simple functional 
form that makes AM/Mo a factor of 7.5 lower than AE/\Eo\ for 
AE/\E \ < 1 and flattens to 1 for AE/\E \ > 1. The AM/Mo 
defined as the number of stars beyond rt is shown as squares. A 
functional fit is shown as a full line. A one-to-one relation is shown 
as dotted line. The results of the simulations are expected to be 
between the grey shaded regions, since the maximum value for 
AM/Mo is 1, corresponding to the total disruption of the cluster 
and the minimum value is 1/8192, when a single star got unbound 
by the perturbation. 



5.2.2 The cluster disruption time 

lOstriker et ail (|l972T ) define the cluster disruption time due 
to many repeated disc shocks (idd) using the impulsive ap- 
proximation by dividing the initial cluster energy Eq by the 
energy injection per unit time: tdd = \Eq\ / (AE / At) . For the 
cluster energy they assume Eo = —0.2 G M c /rh (Note that 
this implies r\ = 0.4, see Section ^. 1.2[) . Combined with equa- 
tion ([8| and the fac t that clusters have tw o disc passages per 
orbital period (P). lOstriker et ah! (|l972l) express tdd as 



3GM C P V 2 



20<? : 



2 ri 



(18) 



This time scale is defined as the time needed to unbind the 
cluster by (periodically) injecting energy in the cluster by 
disc shocks. Note that here we have assumed r h = r 2 . This 
is true only for very concentrated clusters. 

This expression can be rewritten to derive a disrup- 
tion time due to periodic spiral arm perturbations (ids)- We 
explicitly use the mass loss per spiral arm passage (equa- 
tion [T7J, which we approximate as AM/M = f(AE/\E \), 
with / ~ 0.4 (Section 15. 1[1 . The expression for tda can then 
be derived from 



_ At \Ep\ 
ds ~ / AE' 



(19) 



with AE from equation (|13[) (with the adiabatic correction 
factor) and At = tdrift(-R)- Using the definition of Eo from 
Section I4T21 we can write AE/\Eq\ as 
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are shown as a dotted line and as a full line, respectively. 



AE/\E \ = 



2ff m r 2 2r h 
3Fi ax r,GM c 



„3 
' h 



3V^ ax nGM c rl 



(20) 



where we included a 

„3 



Mc/K 
r 



factor 

in equation ([18}. 



" 2 /r" 2 to get a dependence on 
For a cluster with Wo ~ 5 



as 

2 /r 2 ~ 2. 

The adiabatic correction factor (A w (x)) of equation (|20|l 
is calculated using the results of Section [4] We use the pa- 
rameters for the Wo = 5 cluster as described in Sections l4.1.3l 
and 14.21 that is, r\ = 0.41 and u = 0.58. We consider 
clusters with M c = 1O 2 M , 10 3 M and 10 4 M all with 
rh = 3.75 pc. The constant radius is based o n recent obser - 
vations of young clusters in spiral galaxies l|Larsenl l2004h . 
Using the scaling relations of Section 14.1.21 the values of 
to for these clusters are 0.046,0.14 and 0.46 Myr" 1 , respec- 
tively. In panel A of Fig.|8]we show the result for the impul- 
sive approximation, that is A(x) = 1, in the top panel and 
Weinbergs correction factor (A w (x)) in the bottom panel. 

In panel B of Fig. [5] we show the energy gain of equa- 
tion pop using values of A(x) from panel A. For the cluster 
with M c = 10 2 Mq inclusion of the correct A vl (x) does not 
make much difference. For the cluster with M c — 1O 3 M0, 
A w (x) decreases to 0.2 at Rcr and for the most massive 
cluster (M c = 10 4 Mq ) the effect of the spiral arm passage 
is largely damped for almost all values of R. This is because 
t cr is shorter in a more massive cluster, assuming a constant 
radius. 

If we translate AE/\E \ to AM /Mo with equation fTTJ) 
we prevent that the mass loss becomes larger than the clus- 
ter mass. Panel C shows AM /Mo in the impulsive approxi- 
mation (top) and with the adiabatic correction factor (bot- 
tom). The cluster with M c — 10 2 Mq can get completely 
unbound by a single spiral arm passage (AM/Mo=l). 



An expression for t^ B can be obtained by combining 
Eqs.[l9]and[20] 



tds(R) 



3V 2 



1 GM C 



4ffiL f A w (x) 



r 3 
' h 



(21) 



where V ma x ^ iV d rift(-R)| (equation [[} and i drift = 
27ti?/(m| Vdriftl), with m the number of spiral arms in the 
galaxy. In panel D of Fig. [8] we show the result of ids as a 
function of R for a two-armed spiral galaxy (m = 2). For 
clusters inside the corotation radius, i.e. R < Rcr, the value 
of ids is nearly constant. This is because we have assumed the 
arm parameters not to vary with R (Section l2.2p . The vari- 
ables that are a function of R in equation (|2ip are V max and 
tdrift- The product Vm ax tdrift from equation (|21fl varies with 
R approximately as R\l — R/Rcn\, since tdrift oc -R/|Vdrift| 
and Vni ax ~ | Vdrift I (equation [TJ. The product R\ 1 — R/Rcr\ 
increases from at R — to a maximum at 7?cr/2 to de- 
crease again beyond R = J?cr/2. Since at Rcr there is never 
a spiral arm crossing ids is infinite there. 



6 DISCUSSION AND CONCLUSIONS 

Much work has been done on the combined effect of stellar 
evolution, a realistic stellar mass function and a spherically 
symmetric Galactic tidal field to understand the dissolu- 
tion of star clusters. The resul ting cluster disruption times 
of these models (for exampl e Chernof f & Weinberg Il990 



. J90l; 

Takahashi fc Portegies ZwartJ |2000| ; iBaumgardt fc Makinol 
200St ) are longer than the obse rved values for ope n clus- 
ters in the solar neighbour hood (Lamers et al.|[200a ) and in 
the central region of M51 (|Gieles et al.ll2005h . Our study is 
aimed at understanding the effect of perturbation by spiral 
arm passages and to see whether they can contribute to this 
difference. 

From the age distri bution of open clust ers in the solar 
neighbourhood clusters, lLamers et all (|2005h derived a dis- 
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Figure 8. Evolution of the cluster parameters with R in a two armed grand design spiral galaxy. In the top row results are shown 
using the impulsive approximation, that is A(x) = 1, and in the bottom panel we used A vr {x) of equation H15I I. The adiabatic correction 
factor is shown in panel A. Panel B shows AE/\Eq\ for a single shock at different R. The peak around Rcr is because of a minimum 
in Vinax around Rcr- From AE/\Eq\ we derive AM/Mo through equation H17H and the result is shown in panel C. The disruption time 
at different R, t^ s (R) (equation fSTJ , is shown in panel D. The results for clusters with masses of 10 2 M Q 10 3 M and 10 4 M are shown 
with dotted, dashed and dotted-dashed lines, respectively. The hashed region indicated the region where t^ s is smaller than tdrifti an d 
therefore td s from equation l|2f \ can not exist. 



ruption time of 1.3 x 10 3 Myr for a 10 4 Mq cluster, which 
they refer t o as £ 4. A similar analysis was performed by 
iGieles et all {2005) using the age and mass distribution to 
derive the disruption time of clusters in the central region 
of M51. They find U ^ 150 ± 50 Myr for clusters with 
R < 6kpc, that is, within Rcr (Section These ob- 
served disruption time scales can b e compared to the re- 
sults of lBaumgardt fc Makinol (120031 ). They derive the dis- 
ruption time for clusters on a circular motion in the Galaxy 
and consider a realistic stellar mass function and stellar 
evolution and the tidal field of the Galaxy. They used a 
smooth analytical description of the Galaxy potential and 
find £4 = 6.9 x 10 3 Myr for clusters in the solar neighbour- 
hood, that is, almost a factor of five longer than what is 
observed. For the parameters of M51, the predicted value of 
t4 is ten times longer than the observed value. 

Neither the spiral arm perturbations on a simple cluster, 
without mass function and stellar evolution (Section l5.2.2p . 
nor the effect of the Galactic tidal field on a realistic cluster 
can, on their own, explain the observed disruption time. We 
can compare the disruptive effect of the spiral arm pertur- 
bations to the observed values. In Fig. [9] we show the tds 
as a function of cluster mass M c for the solar neighbour- 
hood (top panel) and for a grand design spiral galaxy, such 
as M51 (bottom panel). The parameters of Sections [2] and 
15.2.21 are used, that is, the number of spiral arms is four 
(m = 4) and R ~ 0.9i?cR, with Rcr the corotation radius, 
for the solar neighbourhood and m = 2 and R ~ 0.5i?cR 
for the clusters in M51. The full lines represent the resulting 
disruption times using the impulsive approximation and the 
dashed lines show the result for t ds when the correct adia- 
batic correction is applied. As we showed in Section 15.2.21 



the impulsive result applies to clusters with M c < 3x 10 3 Mq 
in the solar neighbourhood, because of the long crossing time 
t cr of stars in such clusters. For the cluster in a grand design 
spiral at R = 0.5-Rcr the adiabatic correction is important 
for clusters of slightly higher masses, since at that location 
Vdrift is higher (equation [T|. In the solar neighbourhood tds 
becomes constant for M c < 10 2 Mq . This is because a single 
arm crossing can destroy clusters with such low masses and 
tds becomes e qual to Adrift ■ The di srupt ion times from the 
observations of lLamers et al.l l|2005h and lGieles et alj (120051 ) 
for the solar neighbourhood and M51 are indicated as hashed 
regions in the top and bottom panel, respectively. The area 
of the hashed region represents the errors in the observa- 
tions and the mass range of the observed cluster sample. 
Note that both studies found an empirical mass dependence 
of tdisOcMc - 6 . 

For the solar neighbourhood the disruption time due to 
the spiral arm perturbations is about an order of magnitude 
higher than the observed result. For M51 the difference is 
almost three orders of magnitude. This suggest that spiral 
arm perturbations contribute little to the disruption of the 
(low mass) open clusters in the solar neighbourhood and can 
definitely not ex plain the observed short disruption time of 
clusters in M51 l|Gieles et al.l l2005). In Section [5] we argued 
that the contribution of the stellar arms, which we have 
ignored, is even less than that of the gaseous component. 
This because the time it takes the cluster to cross of stellar 
arms is much longer than the crossing time of stars in the 
cluster. This causes the effect of heating by the stellar arms 
to be largely damped adiabatically (Section |3J. 

There should thus be a further disruptive agent, which 
we have not considered so far. The central region of M51 
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Solar neighbourhood: m = 4, <R> = 0.9R R 




Grand design: m = 2, <R> = 0.5R R 




M c [M ] 

Figure 9. Disruption time due to spiral arm passages (t^ s ) for 
different M c in the solar neighbourhood (top panel) and in a 
gra nd design spiral galaxy (bottom). The full lines show the 
result using the impulsive approximation, that is, equation 1210 
with A(x) = 1. The dashed lines show t^ a with A w (a;) (Fig. [SJ - 
The shaded areas are the observed disruption times and mass 
ran ge of the observed clusters from lLamers et all ||2005| ) (top) 
and lGieles et alj ||2005|1 (bottom). 



has a ric h population of giant molec ular clouds (GMCs) (for 
example iGarcia-Burillo et al.lll993h whose velocity with re- 
spect to the cluster is relatively small (~ lOkms - ). Even at 
the relatively large distance from the Galactic centre of the 
solar neighbo urhood there are GM Cs with masses > 10 6 M 
(for example ISolomon et al.|[l98^ . 
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